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ABSTRACT 


A procedure was developed to aid in the testing of a 
two~sided composite null hypothesis about the mean of a 
normally distributed random variable for situations where 
either the population variance is known or unknown. The 
procedure was designed to eliminate the requirement for 
iterative type solution techniques normally used in deter- 
mining the acceptance or rejection region of the subject 
hypothesis. This thesis provides guidelines, curves, and 
tables which will aid in testing a two-sided composite null 
hypothesis. Provisions were also incorporated into the 
procedure to permit testing of hypotheses about the relative 


displacement of the coefficient of variation from zero. 
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SYMBOLOGY 


Probability of a Type I error. 

Probability of a Type II error. 

Critical value defining the rejection or acceptance 
region centered around Uy (o known case) or zero 

(o unknown case). 


The lower bound (Cy) and upper bound (Cp) on the 
critical region. 


; 

xX. 
i=1 + 
—— . 


The true population expected value (mean). 


A specified value of the true population mean. 


Specified lower (d.) and upper (dp) limits for 
the true populatiof mean (uu). 


Population variance (co = population standard 
deviation). 

Sample variance used to estimate a*, specifically; 
2 


a2 n (x; - X) 
ay + dp 
2 . 
sample size. 


c vn 
=. 


§ vn 
a 


Coefficient of variation (77). 


A specified value of the CV. 


The null hypothesis of interest. 


hm > 


DF 


$(x] 


The alternative hypothesis to H,- 
The tail probability that a value equal to or 
More extreme than the observed X would have 
resulted if Hy were in fact true. 
Desired maximum allowable variation of the true 
mean from some specified value (u_) in terms of a 
portion of the population standard deviation 
(i.e., [u-ull < ro). 
The statistic 

Yn (Xx - uy) 


i? 


Degrees of freedom. 


The non-centrality parameter of the non-central 
T distribution, \ = r/n. 


Normal cumulative distribution function. 


INTRODUCTION 


A. BACKGROUND AND PURPOSE 

A procedure was developed to aid in testing a two-sided 
composite null hypothesis about the mean of a normally dis- 
tributed random variable. The intent was to eliminate the 
need for iterative solution techniques when solving for 
critical values manually. In situations where computer 
assistance is inappropriate or unavailable, there is a tendency 
for the analyst to shift from the interval test requirement 
to the computationally easier two-sided simple hypothesis 
(i.e., that the mean is equal to the midrange value of the 
interval). By employing the procedure proposed in this thesis, 
the time required to test a two-sided composite null hypothesis 


by hand can be significantly reduced. 


B. SCOPE 
Procedures were developed for the following test requirements. 


1. For Testing That d,; <u < dp 


For situations where the population variance is known, 
a method was developed for testing that the population mean 


is between two values (dy and dp). 


2. For Testing That {u-u.| < ro 


For situations where the population variance is unknown, 
a method was developed to test that the absolute value of the 


deviation of the population mean from some specified value (uy) 


is less than or equal to a specified portion (r) of the 
population standard deviation. By a slight modification, 
this procedure can be used to test if the coefficient of 
variation is between zero and some specified value. 

Each procedure consists of using relevant situation 
Parameters (i.e., 5, n, 6, a, etc.) to determine the critical 
region using either interpolation curves or tabulated data. 
The procedure was designed to be versatile with respect to 
uSing Type I (a) and Type II (8) error level control, or 


using tail probabilities (P,) for drawing statistical 


a i 
inference. 

To simplify use of the procedure, the paper is 
organized so that the main body contains only 1) minor 
theoretical context, 2) detailed procedural guidelines, and 
3) example applications. The detailed derivation of the pro- 
cedure is provided in Appendix C for readers so inclined. 
Appendices A and B provide curves and tables used to find the 
associated rejection or acceptance regions and the tail area 
probability. The curves are provided to simplify critical 
value interpolation for what is considered to be the majority 
of applicable parameter values. However, the total breadth 
of parameter ranges considered in this paper is covered in 
the tabulated data of Appendix B. It should be noted that 
all explanatory figures in the main body of this paper are 


schematic in nature and do not represent scale drawings. 
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C. SYMBOL AND HYPOTHESIS TESTING CONVENTIONS PECULIAR 
TO THIS THESIS 


Whenever possible, symbol and concept conventions which 
conform to standard hypothesis testing terminology were 
employed. However, because of requirements which are peculiar 
to this procedure, it was convenient in some instances to 
use terminology that may be new to the reader. The following 
symbol and concept conventions apply throughout this Thesis: 

1. The Null Hypotheses (H,) will always be bounded 
intervals such as |i - Mo & 6 or |p ue! < ro. All alterna- 
tive hypotheses are the complements of such intervals. 

2. The level of significance is customarily a pre- 
assigned upper bound on the size of the critical region 
(l.e., that subset of the sample space reserved for rejection 
of H,)- This will hereafter be referred to as a error level 
control. 

3. By reversing the role played by the null and alterna- 
tive hypotheses, it is possible to define an acceptance 
region for the originally formulated interval null hypothesis 
with a specified preassigned upper bound (8) on the size of 


a type II error. This procedure is defined as 8 error level 


control. 


II. RESULTS AND DISCUSSION 


A. CASE OF KNOWN POPULATION VARIANCE 
l. General Remarks 
In order to test the null hypothesis that the popu- 
lation mean falls between two values (i.e., d, S'S dp), 


an iterative solution of equation 1 is required. 


where, 


o is the known population standard deviation, 


Cc is the value of the statistic x which defines 
the critical region centered around 


i qs + de 
Ug is aaa, 
a is the probability of a Type I error, 
da - dy 
56 is equal to ak ace and 


>’ is the integral of the standard normal probability 
density. 


To determine a rejection region for the subject 
hypothesis, an initial value of C would be guessed and the 
right side of equation 1 calculated. A new C would then be 
used iteratively until the desired accuracy in the value of 
a had been obtained. This solution process is not difficult 
or time consuming with the assistance of a computer. However, 
when performed manually, such a process may require more effort 


than is warranted by the scope of the project. 
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The procedure developed in this thesis will provide 
the individual with a method for testing a two-sided com- 
posite null hypothesis which circumvents the previously des- 
cribed iterative solution technique. The procedure derivation 
ia provided in Appendix C. In general, it consists of 


oyn 
normalising equation 1 by reparameteriszing the terma ee 


and on as the parameters "a" and C*, respectively. Appro- 
priate values of C* are then obtained according to the pre- 
scribed procedure which follows, 


2. Procedure Description 


The procedure for testing that ay < ) * de {sa out- 


’ 


Lined in Figure 1 (for a and ¢ error control) and Figure 2 


(for estimating Pip) This procedure was developed for 


situations where asymmetric critical regions are of interest. 
A summary of the procedure follows. 

a. For Controlling The «a Error Level 
wn 
rej 


First calculate a «* . If the value of “a 


ia lesa than or equal to 1.7, the associated value of C* 
is Obtained from Figure A=1l of Appendix A. The lower (cy) 
and upper (Cy) bounds of the rejection region are then 


defined using equation 2. 


If the calculated value of "a" is greater than 
1.7, the relationship between “a" and C* becomes linear and 
equation 3 is used to define the rejection region (Ky 1s 


obtained from Figure 1). If x » Cy or X > Cp then reject 
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a eee 


the hypothesis with a probability equal to gq of rejecting 


when in fact the hypothesis is true. 


(a + Kp) 9 (a + K,)9 
(3) Cy i se eee and Ce = + eee 
n n 


b. For Controlling The 8 Error Level 

First calculate a = a If "a" is less than or 
equal to 3, the associated C* value is obtained from Figure 
Aw-2 (.2 < a < 1.6) or Figure A~3 (1.6 < a < 3) of Appendix A. 
The lower (C)) and upper (Cp) bounds of the acceptance 
region are then defined using equation 2. 

If "a" is greater than 3, the relationship 

between a and C* becomes linear and equation 4 is used to 
define the 8 error level acceptance region. If Cy Ss es Cp: 


the hypothesis that the mean is between day, and dp is accepted 


with a probability of falsely accepting equal to 8. 


(a - K,)9 (a - Kp) 9 
(4) cS # i), «a a oh ee 
L ° “n R ° vn 


c. For Estimating Tail Probability (Pp) 

The exact value of the tail probability (P,) can- 
not be calculated since it depends on the actual value of the 
population mean (uu). However, the procedure depicted in 
Figure 2 will provide an estimate of the maximum possible 


value of P,, for a corresponding value of x equal to or more 


T 
extreme than that observed. The procedure for estimating Pr 


consists of first calculating the values a = beh and 
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My 
ae 


Calculate a= 


Go to Append.A(Fig. 
for .26a<1.6 or Fig.3 
for 1.6¢a<3). | 


Go to Append.A 
Fig. 1. 


Calculate C at desired error Calculate CL & Ce: 


level curve and intersection 


Sik 
of vertical line projected Gingee As Ee 


upward from calculated value ‘a’ Boe mm 


Calculate CL and C 
*« 
For 8 control 


eons 
CL tua 2 3 Catuyt 


For a contro] L 7 


Ca 
vn 


, if XC, or X>C, then If Cs Xs Cy then accept 
| . y 2 
[ren H(d) < us dp) Hod) su < dp) with an 8 
jwith an a probability probability of falsely 
lof falsely rejecting accepting 


SYMBOLS 


d, & d 


=Lower and upper test interval 


ee values. 
§ = -5(dp- qd, ) 
Me™ -5(d) + dp) 
C.F Defines lower bound of critical 
region. 
Cp= Defines upper bound of critical 


region. 


Figure 1: Algorithm For Determining a or £8 Critical Regions 
for Known o 


Yes No No Yes 
error of = .02 
in P, acceptable? f al 
Go to Append. A Calculate Pr as: Calculate Pr from: 
(Fig. 4 for Osas1.5 3 es. ee ee 
or Fig. 5 for Pr 1-aCatC” }tofa-C J, Pre ota-C ] 


where > is the integral 
of the standard normal 
density. 


1.5ea62.3) and deter- 
mine point where 
vertical line from 
‘a' value and 
horizental line 
through C™ intersect. 


Interpolate between 

Sagi to determine 
(Note: these curves 

be plotted log 

and should be inter- 

polated accordinaly). 


Figure 2: Algorithm For Determining The Tail Probability 
For o Known 


(x - ug) Yn 


ce = | ~ . These values are then used to estimate 


Pp by employing the interpolation curves of Appendix A or 


by direct calculation using the equations provided in Figure 2. 


' 


3. Range of Parameter Values Considered 


The procedure itself is suitable for all possible 
parameter values. However, a restriction on the range of 
parameters covered by the curves and tabulated data was 
necessary. The range of parameter values was selected ina 
manner which would provide adequate coverage for the majority 
of conceivable real-world situations. 

a. Curves 

l. For a Error Level Control. Parameter ranges 
covered are a(.01, .025, .05, «075, «1, <5, «2;) and «25) 
and a(0 to 1.7). 

2. For 8 Error Level Control. Parameter ranges 
covered are 8(.0l1, .025, .05, .075, .1, .15, .2, and .25) 
and a(2to 3). 


3. For Estimating Pine Parameter ranges covered 


are a(0 to 2.8) and C*(.04 to 5.2). 
b. Tables 
l. For a Error Level Control. Parameter ranges 
covered are a(.001, .01, .025, .05, .075, .1, .15, .2, and 
-25) and a(0 to 3). 
2. For 8 Error Level Control. Parameter ranges 
covered are §(.001, .01, .025, .05, .075, .1, .15, .2, and 


-25) and a(0 to 5). 


4. Example Applications 


a. Example 1 

In the production of a certain structural element, 
a new size has been established as marketable. The expected 
size of the new element is to be between 16.5 and 17 cm in 
order to satisfy engineering and management requirements. 

For example, if the expected size is below 16.5 cm the 
rejection rate of unacceptable elements will be too high. 
If the expected size exceeds 17 cm, more material is used 
than is necessary. 

For the first 25 units produced, the average 
size is 17.03 cm. It is known from previous experience that 
the element size is normally distributed with a standard 
deviation equal to .2. In addition, no functional relation- 
ship between the mean and standard deviation of the element 
size has been observed in the range of 10 to 20 cm. 

The decision maker is interested in determining 
if the hypothesis that the expected size is between 16.5 and 
17 cm can be rejected with a Type I error equal to .05. The 
test is conducted with 5 equal to .25, Wy equal to 16.75, and 


(.25) (5) 
aa mea 


a equal to or 6.25. Since a is greater than 1.7, 


the rejection region is defined using equation 3, specifically: 


(6.25+1.645)(.-2) 2 16,434 
cide hi ean ae 


Cy = 16.75 = 


and 


Dah) sult ees PlF nA oak ate, acs oh 
erg Rts Cw UE AS rn Oem 5 ase irk eS 


<— _ - 


-m z x a eS ee 8 ee 


Cp = 16.75 + .3158 = 17.066. 
Since x is neither greater than Cp or less than C,, the } 
hypothesis cannot be rejected with a Type I error equal to | 
.05. | 


b. Example 2 


From a sample size of 81 shipping containers 
selected at random, it is necessary to know if the expected 
container weight (which appears to be normally distributed) 
is between 75 and 75.5 kg. The sample average weight and 
variance are 75.22 and 1.44, respectively. A Type II error 
equal to .01 is desired. Since the sample size is relatively 
large, a reasonable estimate of the population variance can 
be obtained using the sample variance. The test is conducted 


with 6 equal to .25, u_ equal to 75.25, and a equal to 


° 
a or 1.875. Since the value of a is less than 3, 

a C* value equal to .0726 is obtained from Figure A-3 of 

Appendix A. Using this value of C*, the values of Ch and 

Ce obtained from equation 4 are 


Cc = [S22 = jess = S00S75 = 15.24 


(.07263) (1.2) 
9 


and 


Cp = 75.25 + .00975 = 75.26. 


Since x is not between Cy and Cre the hypothesis that the mean 


is between 75 and 75.5 cannot be accepted with a Type II 
error equal to .0l1. However, it may be possible to accept 
the hypothesis at a somewhat higher risk (8) of falsely 
accepting. Trying 8 = .05, the new C* obtained from Figure 


A-3 of Appendix A is .35 and accordingly Cy = 75.2 and Ce =) 7563's 


Since x is between Cy and Cpr the hypothesis could be accepted 


with a probability equal to .05 of falsely accepting. 
c. Example 3 
For the situation of Example 1, the decision 
maker wants to know for what values of a the hypothesis 
16.5 < u < 17 could be rejected. The answer can be found by 
estimating the maximum tail probability associated with the 
observed x using ee oe described in figure 2 with 
x-u,)yn 


a = 6.25 and C* = a = 7 Since a+ C* > 3.5; 


the estimate of Po would be P,, = $({.25-7] = .227. Therefore, 


= 
an a probability greater or equal to .23 would be required 


before the hypothesis could be rejected. 


B. CASE OF UNKNOWN POPULATION VARIANCE 

1. General Remarks 

A procedure for situations of unknown population 

variance which would parallel that of the previously dis- 
cussed variance known procedure cannot be developed without 
several stage sampling such as Stein's Procedure (reference l). 
However, the previously developed methods can be extended 
directly to the case of unknown variance provided the null 


hypothesis is scaled by the standard deviation (see below) 


and the noncentral t distribution is used in place of the 
normal. Details follow. 


a. A procedure was developed to test the hypothesis 


lu- HOGI 


|u-u,| < ro or equivalently < r where 


uw is the true population mean, 


u is some specified value of interest near 
which the mean is to be tested, and 


r is the portion of the population o for which 
the test is to determine if the true mean 
is within ro units of Wo: 

b. The procedure can be used to conduct tests about 
the displacement of the coefficient of variation (CV = aT? 
from zero. 

2. Procedure Description 

The procedure for testing that |u- us| < ro is 
provided in figure 3. For testing that CV < CV. figure 4 
should be used. Finally, the method for estimating Pr is 
provided in figure 5. 

As the sample size gets large, the ratio =. converges 
to the test value of r. To take advantage of ile. vekots, 
the curves for a error level control are presented in terms 
of the value &.. instead of just C. This provides an indica- 


n 
tion of the value of r and the sample size n for which the 


3? unknown curves become close to those of the 3? known case. 
However, since the method of estimating Pr requires repeated 
table entry and interpolation, the critical values presented 
in the tables of Appendix B are in terms of C directly. A 


summary of the procedure follows. 


! 


a. For Testing That |. ~ ugl < Fo 


Select the desired a or 8 error level. Consult 
the corresponding curves of Appendix A or tables of Appendix 
B. If the curves are used, multiply the value obtained for 


&. by the sample size (n). If the tables are used, determine 


the corresponding value of C directly for the subject values 
of r and DF. Calculate the statistic 
: 2 ae a 
Vi \ 
If the test is to control the a error level then reject the 
hypothesis, with a Type I error equal to a, if the value of 
r is greater than the C value determined from the a curves 
Or tables. For a test to control the 8 error level accept 
the hypothesis with a Type II error equal to 8 if the value 
of r is less than or equal to the C value obtained from the 
1-8 tables. 
b. For Testing That CV < CV, 
To test the null hypothesis that the coefficient 
of variation does not exceed CV, at the a level of significance, 
One need only apply a role reversal technique to the above 


Procedure. The set equivalences 


wie coy) = eels hy 


(CV <= cv.) = lu ° 


illustrate the fact that the hypothesis of interest is one 
of the alternate hypothesis treated previously, and can be 


Managed by reversing the role and interpretation of a and 


8 error level control from the previous test for 
[u-ul < ro. 

l. For Controlling The a Error Level. Obtain 
the C value corresponding to the subject values of r (equal 
to co and DF from the table of Appendix B associated with 
the oni l-a. Calculate the statistic 


If r < C, reject the hypothesis that CV < CY. with a Type I 
error equal to a. 

2. For Controlling The 8 Error Level. Obtain 
the C value corresponding to the test values of r and DF 
from the curves of Appendix A associated with the value 8 
(i.e., for 8 = .01, consult the curves for an a = .01). 
Calculate the statistic r. If r > C, accept the hypothesis 
with a Type II error = 3. 

c. For Estimating Po 
The actual value of Po associated with the observed 


value of r will depend upon the true value of 


GE usdyn 
i— 


To obtain an estimate of the range of possible values for Pre 
the procedure detailed in figure 5 is used to calculate a 

value for Po corresponding to r set at zero and r set at the 
test value (r). The actual value of Po will be between these 


two extremes. 


Calculate 


For aerror control 


Go to the appropriate figure 
of Append. A for the associated 
a level (Fig. 6-15) or cables 
of Append. 8 for parameter 

values not covered by curves 

of Append. A. Project vertically 
from the r scale to the appro- 
priate OF curve then horizentally 
to the C//n_ scale. Multiply by 

vm to obtain C. 


If r> C, reject the sypothesis 
|Hy with a probability = a of 


jrejecting when in fact it is at 


Figure 3: Algorithm For Testing The Hypothesis A Reb Ss se 


For g error contro] 


Go to the appropriate cable } 
of Append. 8 associated with | 
| the desired s level(i.e. for 
B =.]1,go to the C value table | 
associated with a,Pryles eo9i) 


! 

| 

Project down the correspond ing| 
|r value column and across the | 
| respective DF row to determine) 
| the correct C value (Wote: if | 
| the r value being tested is | 
between two of the table's | 


columns ,interpolation is 
necessary.) 


| 


If r < C accept the hypothesis 
Hy with a probability = 3 of 


accepting when in fact ic is 
\false. 


Calculate r= 


For controlling a error jevel For controlling 8 error level 


Go to the appropriate 3 table 
of Append. 8 (Tables 3 —» 41) 
which correspondsto an3 volue 
equal to the desired a level 
(i.e. for a=.1, go to the C 
critical value cable associated 
with 1-3 =.9). Project down 
the corresponding r value 
column and across the respec- 
tive OF row to determine the 
correct C value. { Note : If 
the r value being tested is 
between two of the vable’s 
columns , interpolate for C by: 


If r < C, reject the »ypothesis 
CV ECV with probability 2 


of rejecting when in fact it is 


High” “Low 


} 


Go to the appropriate a figure 
of Append. A (Fig.6-—15) which 
corresponds to ana value equal 
to the desired 3 value(i.e.for 
8 =.01, go to the a=.0] curves 
of figure 6 or 7). Project 
vertically from the r scale to 
the appropriate OF curve then 
horizontally to the C/»n 
scale. Multiply by Yn to get C 
[ Note : For parameter values 
not covered by the curves ,go 
to the tables of Append. 8 
which provide C_ directly , 
instead of C/vn.] 


If r>C » accept the .ypo- 
thesis CV<CV with 


probability 8 of accepting | 


when in fact it is false. 


Figure 4: Algorithm For Testing The Hypothesis CV < CV, 


25 


than is warranted by the scope of the project. 


i2 


Calculate 


First calculate the estimate of P. for r set equal to Zero as follows: 


Go to Tables 3447.Pick a P_(i.e..5) and go to the associated table 


Project across the respective DF row and down the appropriate 
r column to obtain the correct C value.(Note: If the r beina 
tested falls between two columns,interpolation is necessary. 


Go to a lower Pr value 


P. value table table(i.e. .2). 


t 
(Ties <8). 


= 
oO 
3) 
vr 
“ 
a 
v 
~< 
oO 
nn 


as that © 
been checked on a 
previous step? 
Yes 


[The Pe for this r is between the current table Pr value and the 
previous cable. For a closer approximation interpolate as follows ; 


G 


ae 
(( ee) cP, 


high “low 


Now repeat the 
process to get 
an estimate of 
Pe forr=r. 


Have the ; 
estimates of Pr fo? 


both r=0 and r=r bee 


lies between the two estimates for r=0 and r=r 


Figure 5: Algorithm For Estimating P, For o Unknown 


| 
] 
| 
| 


obtained from Figure 1). If x» GC or x * Cy then reject 
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3. Range of Parameter Values Considered 


Two general levels of r values, along with selected 
incrementation, were chosen to parallel the range considered 
applicable in the majority of situations. For example, in 
most cases where the test of interest is that |p - ugl < ro, 
r would range between zero and l. Conversely, most tests 
concerning the coefficient of variation would require that 
-143 < CV, < 1 (or equivalently, 7 >» r > 1). However, the 
critical values (C) for the total range of r (i.e., o to 7) 
provided in the curves and tables can be used for either 
test. The specific parameter ranges covered were: 

a. Curves 

Curves are provided for a error level control 
for both low and high values of r. 

1. For Low r Values (0 < r < 1). Curves are 
provided £o0r a equal to .01, .025, .05, <1, and <2. 

2. For High © Values (I <x < 7). ‘Curves are 
provided for «a equal to .01, .025, .05, .1, and .2. 

b. Tabulated Data 

The same tables are used to present critical 

values for both a and 8 error level control and for estima- 


ting P To obtain C values for 8 error control, the tables 


r 
associated with a and Po equal to 1-8 are used. The values 
covered by the tables are a, Py, OF 1-8 equal to .0l, .025, 
+OS% wl, why 635 aby «Sy wOy ely 2B) ASy «95, 2975) ahd «99% 


For each value of a, Po and 1-8, three tables are provided 


to cover the range of r from 0 to 7. 


VWwhig store Wt Hei Oe CAALULSEALANY Se A rk 
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4. Example Applications 


a. To Test That |, - Mo < ro 


cn ll i 


It is necessary to determine if the output of a 
new system is between 11.75 and 12 units. The standard 
deviation in the output of the new system is believed (with 
confidence) to be between 2/3 and 3/4 of that of the old 
system (0 = 1.5). Using the conservative estimate of the 
improvement in the system deviation (i.e., 3/4), the test 


(l2 = Bos ae 18), 


parameter r would be set equal to 
For the first 15 units, the values of the sample mean and 
deviation were 11.68 and .3, respectively. The desired level 
of significance (a) is .05. 

The test is conducted with WG equal to 11.875 and 
r equal to 2.517. Using figure A-10 of Appendix A, = is 
determined to be .6, and accordingly, the critical een 
equals (.6)(v15) or 2.324. Since r is greater than C, the 
hypothesis is rejected with a Type I error equal to .05. 

b. To Test That CV < CV, 

A relationship is known to exist between the mean 
and standard deviation of the purity of a certain chemical 
product. A new method of processing is being considered to 
improve the expected purity of the product. The decision maker 
is interested in determining if the same degree of quality 
control would be obtained with the new process as with the 


current method which provide a standard deviation equal to 


-18 of the mean (i.e., a coefficient of variation equal to 


+18). 


for Known oO 


Since a relationship between the mean and standard | 
deviation is known to exist, and only 15 batches of the pro- 
duct have been produced by the new process, the population 


variance is considered unknown. The mean purity of the first 


15 batches was 4.6 parts per million with a sample variance 


equal to .25. The decision maker is concerned that if he 


does not convert to the new method, he will be rejecting, with 

i a probability (a) equal to .05, an opportunity for a better 

product. The rejection region associated with a value of r 

y equal to 5.56 (i.e., art is found by consulting the critical 
value table with 1-8 equal .95. Interpolating between the 
table columns for r equal to 5 and 6, a value of C equal to 
17.69 is obtained. The value of r is calculated to be 35.63. 

§ Since r is greater than C, the hypothesis that the new process 
provides a quality control equal to or greater than the 
current method cannot be rejected with a Type I error equal 


to .05. 


| 
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APPENDIX A 
CURVES FOR DETERMINING aOR 8 CRITICAL VALUES AND TAIL PROBABILITY 
(P 
T 


—_———. 


I. GENERAL 


The curves provided within were developed to simplify 
procedural calculations in lieu of the tabulated data of 
Appendix B. They do not cover the total parameter ranges of 
the data in Appendix B but were designed to maximize interpo- 
lation precision while covering what is felt to be the majority 
of applicable situations. For those cases which are subject 
to parameter values outside the range of the curves, Appendix 
B should be used. Figures A-1 through A-5 are relevant to 
the o? known procedure. Figures A-6 through A-15 deal with 


the o? unknown situation. 
II. ACCURACY 


The accuracy achievable in using the curves vs. using the 
associated tables depends on the interpolation prowess of 
the individual user. However, to provide a general idea of 
the accuracy which may be expected, the following statements 


are considered germane. 


A. < UNKNOWN PROCEDURE CURVES 


a. C* Values For a or 8 Error Level Control 
Accurate interpolation to the third digit is possible 


-005). 


* 
=oeee S interpolated : 


" * 
(1.@., ¢ true 


: , 
REPRE T Amer Caen Tn eter. Tee aay, SS eee en 


b. Po Values 


Because of the semi log scale and interpolation 


between curves, the accuracy in the interpolated Pr values 


is expected to be = + .02 for .1 < P, < .9, and + .005 for 


Pr < .l or Po eee 


B. o? UNKNOWN PROCEDURE CURVES 


a. Curves For The Range of r From 0 to l 


Interpolation accurate to the third digit is 


possible. 


b. Curves For The Range of r From 1 to 7 


Interpolation accurate to the second digit is 


possible. 


ie hice ae Weigh deter Me as 


7) 


C* Values For a Error Level Control With o Known; 
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C* Values For 8 Error Level Control With o Known 
@2< a < 1.63 
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Figure A-4: Tail Probability for o Known with .2< a < 1.5 
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Figure A-5: Tail Probability For o Known With 1.5 Sas 2.8 
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APPENDIX C 


PROCEDURE DERIVATION 
I. KNOWN POPULATION VARIANCE 


A. ERROR LEVEL CONTROL 
The a rejection region associated with the hypothesis 


that the mean is between qs and dp is defined by 


a © Fix <6] Ri > S| 


where the rejection region is chosen to be symmetric around 


Uy (see figure C-l). Since the suprema will occur when 


a = t i 
Uy C=C, Uy-d=d) uy dp=u gts Cer tt 


Figure C-1: Density of x with u at dy ( a =f) 


Lu equals a or dp: the critical region for the specified upper 


bound on the Type I error (a) can be written as 


a = P(x < C, |given ee es P(x > C, | given i = ug7él. 


94 


Since the sample mean (x) will be normally distributed, 


this can be expressed in the form of equation (C-1). 


{C, -(ug- 8) }yn ppt _ tear tig-8) WM 


o Go 


where Zz = tLe {uae$)! ~% N(0,1) 


(C-1) 2 = PL z< 


Substituting C, = Howe and C,, = Hote, equation C-l becomes 


L R 


(C-2) a ae 2 <« LECT 4 > (schon 


or 


a = of — E 0 ] +1 - of — + Cyn 


fo} 


By reparameterizing, the number of parameters necessary 


for solving equation C-2 can be reduced from five to three. 
Letting a = east and C* = SS, equation C-2 becomes 


* * 
(c-3) a= 1+ ofa-C] - ofa +C ] 


Using the result that Cy = Ha > Ss Cp = We tee and 
* 
c= sc", equation C-4 is obtained and is used to define the 
yn 


rejection region. The hypothesis is then rejected, with a 
probability of falsely rejecting equal toa, if the sample 


mean is less than C, or greater than Cp: 
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For each value of a there is a value of "a" such that 

the relationship between "a" and C* in equation C-8 is nearly 
linear (i.e., C* = K, +a). For a > 1.7 the relationship is 
approximately linear (starting with 5 significant digits 

for the worst case with a = .25, and improving as "a" in- 
creases) with slope one for all values of a < .25. The value 
of Kp can be obtained by considering that if ¢{a+cC*] = l, 

then equation C-3 simplifies to equation C-5. This can be 
rewritten to provide K, = @7"[tsal. When this linear relation- 


ship occurs, C, and Ce can be obtained using equation C-6. 


L 


(e-5) «a © ¥e eC 1-7 ce ae tla J 


(Ke+a)yn ; 
(C-6) C ~ ie o and CR > ee o 


B. 8 ERROR LEVEL CONTROL 

In order to determine an acceptance region for the 
hypothesis d, <u S$ dp, with a probability no greater than 
8 of falsely accepting, the alternative hypothesis Ju- ull > 6 
is tested. The a rejection region of the alternative hypothe- 
sis will be the 8 acceptance region of the null hypothesis 
(i.e., d, < u < dy) and will be defined by C, < x < Cy 
(see figure C-2). 

The suprema will occur where u = a, or d,- Thus, the 


upper bound on the Type II error (8) for the hypothesis 


ad, sus d, can be written as 


ews ee 


: 
7 
: 


Rt mc 


Li, R&R 
Figure C-2: Density of x with u at dy ( 8 = shaped area) 


sup = B6(y) = P[C, < xX < C,|given that p = u_- $6) . 
me a ° 
ueH 


This can be expressed in the form of equation C-7. 
(Cp-u.+6 YT (C, -u..+6) 
(C-7) a (H,) = P[z Pes So a - Pz elie, “ib ° %, 


(x -(u, -6) yA 


o 


where Zz = 


i i ay -«- = 
By substituting Ch é Cc and CR Pe and noting 
that z will be normally distributed (i.e.N(0,1)), equation 


C-7 reduces to equations C-8. 


c-a) o(H,) * of EXP). gf UG 
or 
a(H,) = ofa+C)-efa-C] | 


Equation C-8 can be reformulated to resemble equation C-3 
which will permit use of the same solution technique as was 
used for a error control. Specifically, by multiplying 


equation C-8 by -1l, adding 1 to both sides, and setting 


97 


, 
-_ - ———— a ——— end 


l-a (H,) equal to a', equation C-8 reduces to equation C-3. 


The solution of this equation with a' = 1-a (H,) will also 


be the solution for a' = 1-8 (H,)- 


For each value of 8 there is a value of "a" such that 


the relationship between "a" and C* in equation C-8 is nearly 

linear (i.e., C* = K,ta) - For a > 3 the relationship is 
approximately linear (starting with 3 significant digits 

for the worst case of 8 = .0l1, and improving as a increases) 

with slope one for all values of 8 > .01. The value for the 

constant (Ky) can be determined by noting that when 


${a+cC*] = 1, equation C-8 reduces to equation C-9. 


* * -] 
(C-9) a = 1-82 ¢fa-C] or a-C #=¢@ [1 -8 ] 
® h - -)] 
and C sa-Kp where Kp = [) -8 ] 


The value of K, for 8 error control will be the negative of 


E 
the value of Kp for a error control. To allow the same Kp 
value to be used for both a and 8 error control, equation C-10 
would be used for 8 error control when this linear relation- 


ship holds (i.e., a > 3 for this procedure). 


(a- K,) “nm (a-K,) vn 


(C=-10) Cy = Wy ee and CR =u + 


° 10] 


C. ESTIMATING TAIL PROBABILITY 

The actual value of Pr for an observed value of x will 
depend upon the true location of » within the interval da, 
to dp (see figure C-3). Since » is unknown, the exact value 
of Po cannot be calculated. By setting ui = n* 6, the maximum 
possible value of Pr can be Rep ige* from equation C-3 by 
substituting Pr for a, and cet eames for C*. For large 


t values of a+C*, this reduces to equation C-ll. 


(C-11) Pr = of[a-c*] for a+cC* > 3.5 


ite) 
Pr with vau, 


% u -X 
L 0 Uy t(uy=x) 
Pr with ut ule 6 =d) 


Figure C-3: Dependence of Pr on Location of True up Value 


II. UNKNOWN POPULATION VARIANCE 


A. TO TEST THAT |p-u,| < ro 
1. a Error Level Control 
The statistic 
(X- U5) Yn 
Vv oe 


is used to perform this test and has the density of the 


r = 


absolute value of a non-central T, as depicted schematically 
in figure C-4. The associated a error level critical region 
is defined when the non-central location is at r, resulting 
in a non-centrality parameter (\) equal to r/n. The corres- 
ponding critical value (C) is obtained from equation C-12 
and the hypothesis is rejected, with a Type I error equal 


to a, if r is greater than C. 


0 Azrvn C 


Figure C-4: Density of the Non-central T(Absolute Value) 
Distributed # 


(C-12) a = P[ # > C/ given y= rvn] 
2. £8 Error Level Control 


The suprema for the value of 8 will occur when 


} = rvn. Equation C-13 is used to determine the value of C 


which will define the acceptance region with a 8 probability 


of falsely accepting. 
(c-13) PCL | TT) < CC) =. 


By a slight reformulation, the same solution technique 
employed in solving for a critical values can be used here. 
Substituting a' = 1-8 into equation C-12 and solving for C 
will provide the critical value asscciated with the desired 
level of a Type II error = 8. This equivalence is shown in 


equation C-14. 


(C=14) g = PE [Tl] <C) = I+ PC [Tl >C] or PCIT| >C] =1-8 "a 


Because of the equivalence between the C values for 
a and for 1-8, the critical values for a, Pie and 8 can be 
presented in the same table provided that for a desired 
value of a Type II error equal to 8, the table corresponding 


to a = 1-8 is used. 


B. TO TEST THAT CV < CV, 
1. 8 Error Level Control 
To formulate the procedure for testing CV < CV,: 


first consider testing the hypothesis * < r or equivalently 
vou 
° 


| < ©. By setting u, equal to zero, the null hypothesis 
can be rewritten as |f| < xr. Solving equation C-12 with 


We * O andr = |¥“R) wi2 provide the rejection region with a 
6 
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Type I error equal to a. However, rejecting the hypothesis 
that a < xr, when r is greater than C, is equivalent to 
rejecting that CV > i. This in turn is equivalent to 
accepting that CV < CV, with a Type II error equal to the 
value of a. As a result, it is possible to test the hypothe- 
sis that the coefficient of variation is less than some speci- 
fied value (CV) by using the previous test procedure for 
testing that |p - Uo! < ro, provided the following differences 
are observed: 
a. Obtain the critical value (C) from the a tables 
(or curves) associated with a value of a equal to 
the desired value of a Type II error (8) and 


with r equal to ae: 
° 


b. If r is greater than or equal to this value of C, 
accept the hypothesis CV < CV, with a Type II 
error = 8. 
2. a Error Level Control 
The same general formulation philosophy applies for 
a error control of the test CV < Cv, as for 8 error control. 
Again, first consider accepting the hypothesis a < yr with a 


Type II error = 8. This is equivalent to accepting that 


CV > dior Cv). However, accepting the hypothesis CV > 2 
with Type II error = 8 is equivalent to rejecting that 
CV < t (or cv.) with a Type I error equal to the value 6. 


Thus, the desired test that CV < CV, can be conducted for 


controlling the Type I error level by using the test procedure 
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for [u- ug! < ro and observing the following differences: 

a. Obtain the critical value (C) from the 8 tables 
f (i.e., the table associated with the value 1-8), 
} where the value of 8 is set equal to the desired 


F 1 
to Cv" 


| value of a Type I error (a) and r is set equal 
b. If © is less than this value of C, reject the 
hypothesis that CV < CV, with a Type I error 


equal to a. 


B. ESTIMATING TAIL PROBABILITY (Pay) 


The true value of Po will depend on the actual value of 
Seo 
| - =| which, under the null hypothesis, may vary from zero 


to r (see figure C-5). The range of possible values for Pn 
can be estimated by calculating the Pr associated with the 
observed value of r with \ set at the two extremes zero and 


ryn in the solution of equation C-12. 


P(A at we MM 
Pr (a at rvn)= AY 


A 


Rete. BBewne = 
A" MXN NAN IOUS 
420 ABrvn r 


Figure C-5: Dependence of P: on the True Non-centrality Location 


APPENDIX D 


SOLUTION TECHNIQUES AND COMPUTER PROGRAMS 
I. GENERAL 


The same general techniques were used in solving for the 
critical values in both the known and unknown variance 
situations. This consisted of a combined proximity interval 
search and Newton-Raphson technique. A master program was 
used to monitor and control solutions for C values at various 
selected combinations of the respective parameter values. 
Solution accuracy was controlled by stopping the iterative 
procedure when the desired accuracy in the a or 8 error 


level was obtained. The specified accuracy for the case of 


known o* was .000001. In the case of unknown o*, the 


specified accuracy was .00001 for all parameter values 


except the following: 


is @ -Ol1 with 125 degrees of freedom and r = .9 
as. 3 -l with 90 degrees of freedom and r = .9 
3. @ -3 with 125 degrees of freedom and r = 1 
4. a = .3 with 60 degrees of freedom and r = 5 and 6 
5. a= .4 with 50 degrees of freedom and r = 2.5 
The specified accuracy for these latter parameter values 


was .00005. 


II. a? KNOWN SITUATION 


The programs used for this case were SIGMAKWN, CINCR, 
CINCRNEWRAPS, NEWTONRAPS, and NORMPROB. SIGMAKWN was the 
master problem which stepped solutions along the selected 
aor 8 values for each specified "a" value. Controls for 
switching solution subprograms were necessary for parameter 
values whose respective C* were too distant to permit solution 
solely by the Newton-Rapshon method. Equation D-1l was used 


during proximity interval searching (CINCR and CINCRNEWRAPS). 


(D-1) - P(a,a,C*) = kL = @fa + C*] + O[a - C¥] 


Equation D-2 was employed when the Newton-Rapshon method was 


used. 
2.50663 (1-efasc. ,]+fa-C. .]-a} 
‘ -91a P a-l. - 
(D=2) C,*Cy > F C+ ——_— tt ,, 
F/ x % Exp[-.S(a+C,_,) ]+Exp(-.5(a-C;_}) J 


III. 3? UNKNOWN SITUATION 


The programs used for solving for critical values when 
o was unknown were MHUSIGRATIO, MHUSIGRATIOSOLENR, CINCR2, 
CINCRNEWRAP2, NCT, IB, and TRUNERR. MHUSIGRATIO was the 
master program stepping solutions along the selected r values 
for each specified DF. This permitted sole Newton-Raphson 


sOlution once the C for the first r value was determined. 


To solve for critical values associated with testing the 
o* unknown two-sided composite null hypothesis, the equation 
P[r > C] = a (where r has a non-central T distribution) was 
solved for values of r between 0 and 1 using the method 
discussed in reference 2. This method requires use of 


equation D-3. 


ee 


4 (3-1) 
(D-3) P| Tl>c] af feel, y 


1 
(n-1) 
2 


(l-y dy]} 


where x = C+ ( OF + c2) 


To solve the infinite series of the product of a Poisson and 
Beta distribution, the length to which the series must be 
summed in order to maintain a truncation error less than a 
specified amount (c) was determined using the conservative 
estimate provided by equation D-4. Equation D-4 is solved 

by iterating m (the length of the series) until the left hand 


side is less than « (this task is performed by TRUNERR). 


(D-4) Syme 


m+ 


This method of solving the non-central T is quite adequate 


for small degrees of freedom (DF < 30) and non-centrality 


en a ae ee oe ee 


values (\ < 3) but becomes lethargic for nigher values. As 
a result, a second program was written in FORTRAN which used 


the non=central T subroutine of the IMSL packages. 
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